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Abstract 



A comparison of the relative yields of Y resonances in the decay channel in 

PbPb and pp collisions at a centre-of-mass energy per nucleon pair of 2.76 TeV, is 
performed with data collected with the CMS detector at the LHC. Using muons 
of transverse momentum above 4 GeV/ c and pseudorapidity below 2.4, the dou- 
ble ratio of the Y(2S) and Y(3S) excited states to the Y(1S) ground state in PbPb 
and pp collisions, [Y(2S + 3S)/Y(lS)] PbPb / [Y(2S + 3S)/Y(lS)] pp , is found to be 
0.31 lg]g (stat.) ± 0.03 (syst). The probability to obtain the measured value, or lower, 
if the true double ratio is unity, is calculated to be less than 1%. 
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Quantum chromodynamics (QCD) predicts that strongly interacting matter undergoes a phase 
transition to a deconfined state, often referred to as the quark-gluon plasma (QGP), in which 
quarks and gluons are no longer bound within hadrons. Calculations in lattice QCD [1 j indicate 
that the transition should occur at a critical temperature T c ~ 175 MeV, corresponding to an 
energy density e c ~ 1 GeV/ fm 3 . 

If the QGP is formed in heavy-ion collisions, it is expected to screen the confining potential 
of heavy quark-antiquark pairs 0, leading to the melting of charmonium and bottomonium 
states: J/ip, ip(2S), % c , Y(1S), Y(2S), Y(3S), and Xb- The melting temperature depends on 
the binding energy of the quarkonium state. The ground states, ]/xp and Y(1S), are expected 
to dissolve at significantly higher temperatures than the more loosely bound excited states. 
Quenched lattice QCD calculations [3. 4 J originally predicted that the Y states melt at 1.2 T c 
(3S), 1.6 T c (2S), and above 4 T c (IS), while modern spectral-function approaches with complex 
potentials [5] favour somewhat lower dissolution temperatures. This sequential melting pat- 
tern is generally considered a "smoking-gun" signature of the QCD deconfinement transition. 
However, a large fraction of the observed IS yield in elementary collisions is caused by feed- 
down contributions from decays of heavier states (around 50% for the Y(1S) ]6l). Therefore the 
melting of the excited states is expected to result in a significant suppression of the observed 
IS yields, even if the medium is not hot enough to directly dissolve the ground states. 

Observations of ]/ip and ip{2S) suppression between proton-nucleus and heavy-ion collisions 
were reported by the NA38 [7], NA50 d EI and NA60 JTOH fixed-target experiments at the 
Super Proton Synchrotron (SPS), respectively in SU, PbPb and Inln collisions, at centre-of-mass 
energies per nucleon pair (a/snn) of about 20 GeV. The PHENIX experiment, at the Relativistic 
Heavy Ion Collider (RHIC), extended the J/i/> suppression measurements to AuAu collisions 
at y'SNN = 200 GeV iTlTI . At RHIC, bottomonia production becomes measurable [12], though 
with limited integrated luminosities. PHENIX observed that the dimuon yield in the Y mass 
region for minimum bias AuAu collisions is less than 64%, at the 90% confidence level, of the 
value expected by extrapolating the pp yields IIT31 . 

A new era of detailed studies of the bottomonium family in heavy-ion collisions has started at 
the Large Hadron Collider (LHC). The measurement reported in this Letter is performed with 
the data recorded by the Compact Muon Solenoid (CMS) experiment during the first PbPb LHC 
run, at the end of 2010, and during the pp run of March 2011, both at i/snn = 2.76 TeV The 
integrated luminosity used in this analysis corresponds to 7.28 f/b -1 for PbPb and 225 nb _1 
for pp collisions, the latter corresponding approximately to the equivalent nucleon-nucleon 
luminosity of the PbPb run. The excellent momentum resolution of the CMS detector results 
in well-resolved Y peaks in the dimuon mass spectrum. The CMS collaboration has previously 
studied Y production in pp data at \/s = 7 TeV [14J, using techniques to extract the Y yields 
that are very similar to the ones used in the study reported in this Letter. 

A detailed description of the CMS detector can be found in |T5| . Its central feature is a su- 
perconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. Within 
the field volume are the silicon pixel and strip tracker, the crystal electromagnetic calorimeter, 
and the brass /scintillator hadron calorimeter. Muons are measured in gas-ionisation detectors 
embedded in the steel return yoke. In addition, CMS has extensive forward calorimetry, in 
particular two steel/ quartz-fiber Cerenkov hadron forward (HF) calorimeters, which cover the 
pseudorapidity range 2.9 < \rj\ < 5.2. 

In this analysis, Y mesons are identified through their dimuon decay. The silicon pixel and 
strip tracker measures charged-particle trajectories in the range \rj\ < 2.5. The tracker consists 
of 66M pixel and 10M strip detector channels, providing a vertex resolution of ~ 15 f/m in the 
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transverse plane. Muons are detected in the \t]\ < 2.4 range, with detection planes based on 
three technologies: drift tubes, cathode strip chambers, and resistive plate chambers. Due to 
the strong magnetic field and the fine granularity of the silicon tracker, the muon transverse 
momentum measurement (pj) based on information from the silicon tracker alone has a reso- 
lution between 1 and 2% for a typical muon in this analysis. 

In both the PbPb and pp runs, the events are selected by the CMS two-level trigger. At the first, 
hardware level, two independent muon candidates are required in the muon detectors. No 
selection is made on momentum or pseudorapidity, but in the pp case more stringent quality 
requirements are imposed for each muon in order to reduce the higher trigger rate. In both 
cases, the software-based higher-level trigger accepts the lower-level decision without applying 
further criteria. From reconstructed JPsi — > decays, the single-muon trigger efficiencies are 
measured and found to be consistent between the PbPb, (96.1 ± 1.0)%, and the pp, (95.5 ± 
0.6)%, data sets, for muons with pj > 4 GeV/ c. 

In the PbPb data, events are preselected offline if they contain a reconstructed primary vertex 
made of at least two tracks, and a coincidence in both HF calorimeters of energy deposits in 
at least three towers of 3 GeV each. These criteria reduce contributions from single-beam in- 
teractions (e.g. beam-gas and beam-halo collisions with the beam pipe), ultra-peripheral elec- 
tromagnetic collisions, and cosmic-ray muons. A small fraction of the most peripheral PbPb 
collisions is not selected by these requirements, which accept (97 ± 3)% of the hadronic inelastic 
cross section [16 1. For the pp run, a similar event filter is applied, relaxing the HF coincidence 
to one tower in each HF, with at least 3 GeV deposited. This filter removes only 1% of the pp 
events satisfying the dimuon trigger. 

The muon offline reconstruction is seeded with ~ 99% efficiency by tracks in the muon detec- 
tors, called standalone muons. These tracks are then matched to tracks reconstructed in the 
silicon tracker by means of an algorithm optimised for the heavy-ion environment fl7llT8| . For 
muons from Y decays the tracking efficiency is ~ 85%. This efficiency is lower than in pp, as 
in PbPb the track reconstruction is seeded by a greater number of pixel hits to reduce the large 
number of random combinations arising from the high multiplicity of each event. Combined 
fits of the muon and tracker tracks are used to obtain the results presented in this Letter. The 
heavy-ion dedicated reconstruction algorithm is applied to the pp data in order to avoid po- 
tential biases, arising from different tracking efficiencies of the two reconstruction algorithms, 
when comparing the two data sets. 

Identical very loose selection criteria are applied to the muons in the pp and PbPb data. The 
transverse (longitudinal) distance from the event vertex is required to be less than 3 (15) cm. 
Tracks are only kept if they have 11 or more hits in the silicon tracker and the x 2 P er degree of 
freedom of the combined (tracker) track fit is lower than 20 (4). The two muon trajectories are 
fit with a common vertex constraint, and events are retained if the fit x 2 probability is larger 
than 1%. This removes background arising primarily from displaced heavy quark semilep- 
tonic decays. As determined from Monte Carlo simulation of the Y(1S) signal, these selection 
criteria are found to reduce the efficiency by 3.9%, consistent with the signal loss observed in 
both pp and PbPb data. The available event sample limits to 20 GeV/ c the dimuon transverse 
momentum range probed in this study. 

In order to further reduce the background in the Y mass region, only muons with a transverse 
momentum (pt) higher than 4 GeV/c are considered, resulting in a Y acceptance of approxi- 
mately 25% for the \y Y \ < 2.4 rapidity range. This requirement improves the significance of 
the Y(1S) signal in PbPb data and is applied to both data sets. The acceptance of a Y state 
depends on its mass, since the excited states give rise to higher-momenta muons. In con- 
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sequence, requiring higher p h T increases the acceptance for the excited states relative to the 
ground state. In the corresponding analysis performed with the higher-statistics (3.1 pb _1 ) 
7 TeV data JlH, looser criteria were applied {p\ > 3.5 GeV/c and \r/ }l | < 1.6, or p f j >2.5 GeV/c 
and 1.6 < \r/V\ < 2.4), where tjV is the muon pseudorapidity. The stricter {pt > 4 GeV/c) 
requirements used here enhance the Y(2S + 3S)/Y(1S) yield ratio by ~ 60% in the pp data 
at 2.76 TeV. It was checked that, applying the same reconstruction algorithm and the same pj 
requirements, the Y(2S + 3S)/Y(1S) yield ratio is consistent between the 2.76 and 7 TeV pp 
data sets. 

The dimuon invariant mass spectra with the selection criteria applied are shown in Fig. [I] for 
the pp and PbPb data sets. Within the 7-14 GeV/ c 2 mass range, there are 561 (628) opposite- 
sign muon pairs in the pp (PbPb) data set. The three Y peaks are clearly observed in the pp case, 
but the Y(2S) and Y(3S) are barely visible over the residual background in PbPb collisions. 

An extended unbinned maximum likelihood fit to the two invariant mass distributions of Fig. [I] 
is performed to extract the yields, following the method described in [14]. The measured mass 
lineshape of each Y state is parameterised by a "Crystal Ball" (CB) function, i.e. a Gaussian re- 
solution function with the low-side tail replaced by a power law describing final state radiation 
(FSR). Since the three Y resonances partially overlap in the measured dimuon mass, they are 
fit simultaneously. Therefore, the probability distribution function (PDF) describing the signal 
consists of three CB functions. In addition to the three Y(nS) yields, the Y(1S) mass is the only 
parameter left free, to accommodate a possible bias in the momentum scale calibration. The 
mass differences between the states are fixed to their world average values [19] and the mass 
resolution is forced to scale with the resonance mass. The Y(1S) resolution is fixed to the value 
estimated in the simulation, 92 MeV/ c 2 , which is compatible with the resolution obtained from 
both the PbPb and pp data. The low-side tail parameters are also fixed to the values obtained 
via simulation. Finally, a second-order polynomial is chosen to describe the background in the 
7-14 GeV/ c 2 mass-fit range. 

The quality of the unbinned fit is checked a posteriori by comparing the obtained lineshapes to 
the binned data of Fig.Jl] The x 2 probabilities are 74% and 77%, respectively for pp and PbPb. 

The ratios of the observed yields of the Y(2S) and Y(3S) excited states to the Y(1S) ground 
state in the pp and PbPb data are: 

Y(2S + 3S)/Y(lS)| pp = 0.78+°^ ±0.02, (1) 
Y(2S + 3S)/Y(lS)| PbPb = 0.24^ ±0.02, (2) 

where the first uncertainty is statistical and the second is systematic. 

The systematic uncertainties are computed by varying the lineshape in the following ways: (1) 
the CB-tail parameters are varied randomly according to their covariance matrix and within 
conservative values covering imperfect knowledge of the amount of detector material and FSR 
in the underlying process; (2) the resolution is varied by ±5 MeV/ c 2 , which is a conservative 
variation given the current understanding of the detector performance and reasonable changes 
that can be anticipated in the Y-resonance kinematics between pp and PbPb data; (3) the back- 
ground shape is changed from quadratic to linear while the mass range of the fit is varied from 
6-15 to 8-12 GeV/ c 2 ; the observed root-mean-square of the results is taken as the systematic 
uncertainty. The quadrature sum of these three systematic uncertainties gives a relative uncer- 
tainty on the ratio of 10% (3%) on the PbPb (pp) data. 

The ratio of the Y(2S + 3S)/Y(1S) ratios in PbPb and pp benefits from an almost complete 
cancellation of possible acceptance and /or efficiency differences among the reconstructed re- 
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Figure 1: Dimuon invariant-mass distributions from the pp (a) and PbPb (b) data at ^Snn = 
2.76 TeV. The same reconstruction algorithm and analysis criteria are applied to both data sets, 
including a transverse momentum requirement on single muons of pj > 4 GeV/ c. The solid 
lines show the result of the fit described in the text. 
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sonances. A simultaneous fit to the pp and PbPb mass spectra gives the double ratio 

Y(2S + 3S)/Y(lS)|p b p b +0 .i9 / , . \ _i_ n no fsvst \ 

Y(2S + 3S)/Y(lS)|pp " °' 3 -°' 15 (S atj U ' U3 (SyS } ' (3) 

where the systematic uncertainty (9%) arises from varying the lineshape as described above in 
the simultaneous fit, thus taking into account partial cancellations of systematic effects. 

The single muon lower momentum requirement is a posteriori varied from 3 to 5 GeV/ c in steps 
of 500 MeV/ c, and it is found that pj requirements other than 4 GeV/ c provide lower double 
ratios. Fitting the pp and PbPb spectra with free and independent mass resolution parameters 
leads to an increase of the double ratio by 15%. 

To evaluate possible imperfect cancellations of acceptance and efficiency effects in the double 
ratio, a full GEANT4 [20J detector simulation is performed. The effect of the higher PbPb un- 
derlying event activity is considered by embedding, at the level of detector signals, Y(1S) and 
Y(2S) decays simulated by pythia 6.424 [21] in PbPb events simulated with hydjet J22J. Track 
characteristics, such as the number of hits and the x 2 of the track fit, have similar distributions 
in data and simulation. As mentioned above, the trigger efficiency is evaluated with data, by 
using single-muon-triggered data events, and reconstructing ]/tp signal with and without the 
dimuon trigger requirement. The same exercise is carried out with the simulation and it agrees 
with the efficiency measured in data at the 2% level. The track efficiency in the silicon detector 
is measured with standalone muons, applying all selection criteria. The efficiencies in data and 
simulation agree within the 4% statistical uncertainty of the efficiency determined from data. 

The difference in reconstruction and selection efficiencies between the Y states is less than 5% 
and the variation with charged particle multiplicity is less than 10% from pp to central PbPb 
collisions, producing a maximum change of 0.5% on the double ratio. The good agreement 
between single-muon trigger efficiencies extracted from data for the pp and PbPb trigger re- 
quirements, applied to the Y(1S) and Y(2S) trigger efficiencies derived from simulation, leads 
to a negligible effect on the double ratio. The single-muon trigger efficiencies extracted from 
data agree within 1.5% for the pp and PbPb trigger requirements, and the Y(1S) and Y(2S) 
trigger efficiencies agree within 3%, according to simulation: the potential trigger bias on the 
double ratio is negligible. The magnitudes of the statistical and systematic uncertainties on the 
double ratio, respectively 55% and 9%, are significantly larger than the systematic uncertainties 
associated with possible imperfect cancellation of acceptance and efficiency effects. Therefore 
no additional uncertainty from these sources is applied. 

Finally, using an ensemble of one million pseudo-experiments, generated with the signal line- 
shape obtained from the pp data (Fig.jl^), the background lineshapes from both data sets, and a 
double ratio (Eq. [3]) equal to unity within uncertainties, the probability of finding the measured 
value of 0.31 or below is estimated to be 0.9%. In other words, in the absence of a suppression 
due to physics mechanisms, the probability of a downward departure of the ratio from unity of 
this significance or greater is 0.9%, i.e. that corresponding to 2.4 sigma in a one-tailed integral 
of a Gaussian distribution. 

Other studies from the CMS experiment show that the Y(1S) itself is suppressed by about 
40% [23 [ in minimum bias PbPb collisions at a/snn = 2.76 TeV. Since a large fraction of the 
Y(1S) yield arises from decays of heavier bottomonium states [6j, this Y(1S) suppression could 
be indirectly caused by the suppression of the excited states reported in this Letter. 

Production yields of quarkonium states can also be modified, from pp to PbPb collisions, in 
the absence of QGP formation, by cold nuclear matter effects [24J. However, such effects 
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should have a small impact on the Y double ratio reported here. The nuclear modifications 
of the parton distribution functions (shadowing) should have an equivalent effect on the three 
Y states, because their production involves very similar partons, cancelling in the ratio, at least 
to first order. The same should happen to any other initial-state nuclear effect. In principle, 
the larger and more loosely bound excited quarkonium states are more likely to be broken up 
by final-state interactions while traversing the nuclear matter, something extensively studied 
in the context of charmonium suppression at lower energies 11251 . This "nuclear absorption" 
becomes weaker with increasing energy, and should be negligible at the LHC. At RHIC en- 
ergies, the STAR experiment [26] has reported a Y(1S + 2S + 3S) yield in dAu collisions of 
0.78 ± 0.28 ± 0.20 times the yield expected by scaling pp collisions, compatible with the ab- 
sence of absorption. Furthermore, the double ratio presented here would only be sensitive to 
a difference between the nuclear dependencies of the three states and already at much lower 
energies the Fermilab E772 experiment observed ||27|, in proton-nucleus collisions, no such 
difference, within uncertainties, between the Y(1S) and the sum Y(2S + 3S). 

Future high-statistics heavy-ion and proton-nucleus runs at the LHC will provide further quarko- 
nia measurements, which should help disentangle nuclear from medium effects and aid the 
interpretation of the result reported in this Letter. 

In summary, a comparison of the relative yields of Y resonances has been performed in PbPb 
and pp collisions at the same centre-of-mass energy per nucleon pair of 2.76 TeV. The double 
ratio of the Y(2S) and Y(3S) excited states to the Y(1S) ground state in PbPb and pp collisions, 
[Y(2S + 3S)/Y(lS)] PbPb /[Y(2S + 3S)/Y(lS)] pp , is found to be 0.31 (stat.) ± 0.03 (syst), 
for muons of pj > 4 GeV/c and \r/\ < 2.4. The probability to obtain the measured value, or 
lower, if the true double ratio is unity, has been calculated to be less than 1%. 
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